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ABSTRACT 

This paper investigates complex effects in discrete time prey-predator model with harvesting on prey. The fixed 
points and their local stability were analyzed. Harvesting is strong impact on the dynamics evaluation of 
population .To a certain extent it can control the long term stationary density of population efficiently. 
Bifurcation diagram has been obtained for selected range of different parameters. As some parameters varied 
the model exhibited chaos as long time behavior. 
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I. INTRODUCTION 

Among all mathematical models, prey -predator model have received much attention during the last 
few decades due to its wide range of application. The dynamic relationship between prey-predator and their 
have long had been one of the dominant themes due to its universal prevalence and importance. The prey- 
predator interaction has been described firstly by two pioneers Lotka [1] and Volterra [2] in two independent 
works. After them, more realistic prey-predator models were introduced by Holling [3] suggesting three kinds 
of functional responses for different species to model the phenomena of predation. As is well known in the 
theory of population dynamical models there are two kinds of mathematical models: The continuous -time 
models described by differential equations, and the discrete time models described by difference equations, In 
recent years more and more attention is being paid to discrete time population models. The reasons are as 
follows first the discrete time models are more appropriate than the continuous time models when populations 
have non-overlapping generations, or the number of population is small. Second we can get more accurate 
numerical simulations results from discrete time models. Moreover the numerical simulations of continuous 
time models are obtained by discretetizing models. At last discrete time models are rich dynamical behaviours. 
Never less, some of works by {4],[5],[6],[7][8]. Showed that for discrete time prey-predator models the 
dynamics can produce a much richer set of patterns than those observed in continuous time models. [9] analyzed 
prey-predator model incorporating a prey refuge and independent harvesting on either species. Several authors 
improved by the basic Lotka- Volterra model by introducing Allee effect and functional response [10], [11], [12] 
Harvesting has a strong impact on the dynamic evaluation of a population subject to it depending on the native 
of the harvesting strategy the long run stationary density of population may be significantly smaller than long 
run stationary density of population in the absence of harvesting. Therefore while a population can in the 
absence of harvesting to free extinction risk, harvesting can read to incorporation of a positive. 

II. MATHEMATICAL MODEL 

This paper consider the following system difference equation which describes interaction between 
two species prey-predator 

x(n + l) = rx(n)(l- x(n) - ax(n)y(n) - hx(ri) 

y(n + \) = {\-c)y(ri) + bx(n)y(ri) (1) 
Where r represents the natural growth of the prey and a represents' effect of predation on the prey C 
represents natural death rate of predator and b represents the efficiency and propagation rate of the 
predator in the presents of the prey and h represents harvesting effect 
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III. FIXED POINT AND LOCAL STABILITY 

We now study the existence of fixed point of the system (1) particularly we are interested non-negative 
or interior fixed point. To begin with we list all possible fixed points. 

By simple computation of the above algebraic system it was found that there are three non-negative fixed points 

i) E Q (0, 0) is the origin. 

r _(!+&) 

ii) E. ( , 0) is axial fixed point in the absent of predator exist for r>\ + h 

r 

r , * »x . . ■ ■ c , . . * c . r(b-c)-b(l+h) 

ill) £, 2 (X ,y ) is the interior fixed point where x =— y = 

b ab 



IV. DYNAMIC BEHAVIOR OF THE MODEL 

Now we investigate the local behavior of the model (1) around each of the above fixed points. The 
local stability analysis of the model (1) can be studied by computing the variation matrix to each fixed point. 
The jacobian matrix of the equation (1) at the state variable is given by 



J(x, y) = 



r(l — 2x) — ay—h 
by 



-ax 



bx+(l-c) 



(2) 



The characteristic equation of the jacobian matrix can be written as 
A 2 -TrA + Det = 0 



(3) 



Where Tr is the trace and Det is the determinant of the Jacobian matrix J(x, y) which is defined as 
Tr = [r(l - 2x) - (ay + h) + bx + (1 - c)] and Det = [(r(l - 2x) - h)(bx + (1 - c) - ay + aye] 



Hence the model (1) is dissipative dynamical system is if I [(r(l — 2x) — h)(bx + (1 — c)] — ay + aye l< 1 
Conservative dynamical one if and only if 

I [(r(l — 2x) — h)(bx + (1 — c)] — ay + aye 1= 1 and is un dissipated dynamical system otherwise. 

In order to study the stability of fixed points of the model, we first give following lemma which can be easily 
proved by relation between roots and co-efficient of quadratic equation (3) 

Non-linear systems are much harder to analyze since in most cases they do not possesses quantitative 
solution even when explicit solution are available they are often too complicated to provide much insight. One 
of the most useful techniques for analyzing non-linear system quantitavely is the linearised stability technique 
the stability of the system is investigated by obtaining eignvalues of the jacobian matrix is associated with 
fixed points in order to study the stability of the fixed point model we first give the following theorem 

Theorem: Let p(X) = X 3 +BX 2 +CX + D be the roots of p(X) = 0 .Then the following 
statements are true 



a) If every root of the equation has absolute value less than one, then the fixed point of the System is 
locally asymptotically stable and fixed point is called a sink 

b) If at least one of the roots of equation has absolute value greater than one then the fixed point of the 
system is unstable and fixed point is called saddle 

c) If every root of the equation has absolute value greater than then the system is a source. 

d) The fixed point of the system is called hyperbolic if no root of the equation has absolute value equal 
to one .If there exists a root of equation with absolute value equal to one then the fixed point is called 
Non-hyperbolic. 

Proposition: 1 The fixed point E 0 of the system is locally asymptotically stable if r — h < 1 and 0 < C < 2 
;unstable if r — h < 1 and C > 2 otherwise. 
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Proof: By linearizing system (1) at E Q , one obtain the Jacobian 
r-h 0 
0 1-c 



The eigenvalues of the Jacobian matrix at E 0 are \ = r — h and X 2 = 1 — C 

Thus it clear that E Q is locally asymptotically stable when r — h < 1 and 0 < C < 2 , unstable when r — h< 1 
and C> 2 

Proposition: 2 The fixed point E l of the system is locally asymptotically stable if l<r — /?<3and 



otherwise unstable fixed point. 

r-d+h) 

Proof: one can see that the Jacobian matrix for E. 



2 -(r-h) 
0 



a(L + h — r) 



r _ (1 + /j) 

The eigenvalues of the Jacobian matrix at E l are A, =2 — (r — /l) and =^( ) + (l — c) 



Hence £?] is locally asymptotically stable whenl<r — h<3 and br < 



cr 



r-(l + h) 



and unstable when 



cr 



\<r — h<3 and br > 

r-{\ + h) 

Proposition: 3 

The fixed point E 2 (x ,y ) of the system is locally asymptotically stable if 

b(c + h - 4) b 
(c + h)(b-(2+c + h)) r b-c(l+h) 

Proof: 

We now investigate the local stability and bifurcation of interior point E 2 . The Jacobian matrix (1) at E 2 has the 
form 



J(E 2 ) 



1 



cr 



-ac 



b b 
r(b-c)-b(l + h) 



1 



trJ(E 2 ) = 2- 



rc 



(4) 



DetJ(E 2 ) = l-c(X+h) + cr(l - 



(5) 



By solving the characteristic equation (3) at the interior fixed point E 2 the roots (Eigen Values at E 2 ) will be 
_ Tr(J(Ei)) ± ^(Tr(J(Ei)2 - ADetJ{Ei)) 



A L1 =(l-j)±J(l-^f-c((r- T (l + c)))-(Uh) 



rc . 
2b' 



r 

— i 
b 
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Where Tr(J(Ei)) and DetJ(Ei) given by equations (4) and (5). The Eigen Values in numerical can be 
used to classify different type of bifurcation. 

V. NUMERICAL SIMULATION 

In this section we give the numerical simulations to verify our theoretical results proved in the 
previous section by using MATLAB programming. We also confirm the results by visual representation of the 
system for some values of parameters. We provide some numerical evidence for quantitative dynamic behavior 
of the map. Bifurcation diagram are consider in following cases 
[1] Fixing parameters a=2, b=3, c=0.8, h=0 and varying r 
[2] Fixing parameters a=1.5, b=3. c=0.9, h=0.02 and varying r 
[3] Fixing parameters a=1.3, b=3.34, c=0.9, r=2 to 3.6 and varying h 
[4] Fixing parameters a=1.3, b=3.34, c=0.9, r=2 to 3.8 and varying h 

The case (1) bifurcation diagram of the map(l) in x-y plane is showing the dynamical behavior of the prey - 
predator model system as a varying and fixing parameters a=2,c= a=2, b=3, c=0.8, h=0. From Figure (1) one 

can see that the orbit with initial condition (0.4, 0.3) approaches a stable fixed point E 2 for r < 2.69 and a 

discrete Hopf bifurcation occurs at r=2.69.As r increases the interior fixed point becomes unstable through a 
discrete Hopf bifurcation and the behavior of prey-predator model becomes chaotic. It means that for large 
values of natural growth rate of the prey species r, a system converge always to complex dynamics .In the case 
of r > 2 one observes a discrete Hopf bifurcation occurs and complex dynamics behavior begin to appear for 
2<r<3.8. 
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Bifurcation diagram prey population 
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FIG (i) 



Bifurcation diagram prey population 
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FIG (ii) 



Bifurcation diagram prey populat 
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FIG (iii) 



A bifurcation diagram with respect to r in Fig2 while others parameters are fixed as follows a=l,5,b=3 
,c=0.9,h=0.02 .also bifurcation diagram with respect to a plotted in figure 3 others parameters are fixed as 
follows a=1.3, b=3.34,c=0.9, r=2 to 3.6 from Figures (1-3) the change of parameter a affect of stability of 
the system under goes Hopf bifurcation. 
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Bifurcation diagram prey population 
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Bifurcation diagram prey population 



2 2.2 2.4 




FIG (iv) 



FIG (v) 



Bifurcation diagram prey populatio 
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FIG (vi) 



VI. CONCLUSION: 

In this paper we analyzed the complex effects in discrete time prey-predator system. The stability of fixed 
points, bifurcation and chaotic behavior are investigated in this system. The influence of the main parameters ( 

r, a,b, C, K) on the local stability is studied. We observed that under certain parametric condition the interior fixed point 

enters discrete Hopf bifurcation phenomenon. It also been observed that when harvesting activity is not considered then the 
equilibrium of prey-predator model (1) is stable and when harvesting activity of prey is taken into consideration then the 
population size of the predator decreased and the natural stable equilibrium of the model becomes unstable. 
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